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The opening of a gap in single-layer graphene is often ascribed to the breaking of the equivalence 
between the two carbon sublattices. We show by angle-resolved photoemission spectroscopy that 
Ir- and Na-modified graphene grown on the Ir(lll) surface presents a very large unconventional 
gap that can be described in terms of a phenomenological massless Dirac model. We discuss the 
consequences and differences of this model in comparison of the standard massive gap model, and we 
investigate the conditions under which such anomalous gap can arise from a spontaneous symmetry 
breaking. 
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I. INTRODUCTION 

The isolation of single-layer and few-layer graphene 
has triggered a huge burst of interest, mainly motivated 
by the observation of unconventional electronic proper¬ 
ties, which stem from the Dirac-like low-energy electronic 
structure of graphene characterized by a gapless coni¬ 
cal dispersion [lj]. The huge electronic mobility of free¬ 
standing graphene originates from its chiral properties, 
tightly linked with the lack of electron backscattering 
phenomena near the Fermi level [l|. A drawback of this 
characteristic band structure is the absence of an energy 
gap between the Dirac cones, which would be highly de¬ 
sirable for the exploitation of graphene in device appli¬ 
cations. So far, however, the effective employment of 
graphene-based materials in low-energy electronics has 
been hindered by the difficulty of opening a bandgap 
without affecting the electronic mobility. Understanding 
the fundamental mechanisms responsible of gap open¬ 
ing in graphene is thus of the highest importance in the 
perspective of engineering new efficient switch on/off de¬ 
vices. 

From the theoretical point of view, the simplest way 
to open a gap in the conical Dirac-like dispersion of a 
two-dimensional honeycomb material is to induce an in¬ 
equivalence between the two carbon sublattices A and B. 
This corresponds to include a oc a z term in the Dirac-like 
Hamiltonian: 

H = E D i + ^2 hv [k x a x + ky&y] + y & z , (1) 

k 

where v is the Dirac velocity, k is the momentum relative 
to the K point, and where cb are 2x2 Pauli matrices de¬ 
fined in the subspace of the two carbon orbitals for unit 
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cell. The term Ed represents an energy off-shift of the 
Dirac point due to finite doping. The total dispersion re¬ 
sults thus: Ek t ± = E D ±^/(hvk) 2 + A 2 /4, where k = |k|. 
Note that this model predicts for small k a parabolic be¬ 
havior Ek,± ~ -Ed ± (A/2 + h 2 k 2 /2m e s) for both con¬ 
duction and valence bands with an effective mass m e s 
proportional to the band gap A: m e s = A 2 /2v. In order 
to highlight the strict relation between the gap opening 
and the onset of an effective mass, we define this sce¬ 
nario as “massive” (ms) gap model, in contrast with a 
“massless” (ml) gap model that we will discuss below. It 
should be worth to note, in addition, that, in the massive 
gap model, the band dispersion recovers a standard lin¬ 
ear behavior (Ek,± ~ Ed ± Hvk) for Hvk A/2, so that 
the extrapolation of the upper band dispersion overlaps 
the lower band and vice versa. 

Angle-resolved photoemission spectroscopy (ARPES) 
is one of the most direct experimental techniques that al¬ 
low to investigate the energy-momentum dependence of 
the electronic states in solids. Although ARPES can be 
applied also to suspended graphene, as recently demon¬ 
strated in Ref. it is more often applied to study 
the electronic structure of supported graphene. When 
graphene is in contact with a solid material finite dop¬ 
ing effects, as well as modifications (screening) of the 
many-body interactions, are observed. Gap opening is 
frequently reported, with magnitudes that depend on the 
growth process and on the substrate [3Ml7l| . These ob¬ 
servations are commonly interpreted in term of the mas¬ 
sive model discussed above. At a closer look, however, 
the actual evidence of such physical behavior in most of 
the cases is not assessed. A first issue regards the effec¬ 
tive opening of a gap in some epitaxially grown graphene 
samples, where a spectral anomaly at the Dirac point 
was shown to be associated with the signature of plas- 
maronic subbands rather than with a simple gap open¬ 
ing [l§ - l23) . On the other hand, ARPES data for other 
supported graphene layers do not show the characteris¬ 
tic “diamond-shaped” dispersion of the plasnraron at the 
Dirac point fglll4lll7l]. Even in these cases, however, the 
actual experimental band dispersion presents many in- 
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consistencies with the simple massive Dirac gap model, 
pointing thus towards different mechanism of gap open¬ 
ing- 

The ARPES data of Ref. |3. 1{, for instance, reported a 
gap A = 0.26 eV which was initially discussed in terms of 
the massive gap model, and critically revised in Ref. 24]. 
Two main unconventional features were there pointed out 
in regards to the ARPES data: (i) in spite of a large gap 
opening, the conduction and valence bands in Ref. [3||4( 
retained a conical shape, in contrast with the expected 
parabolic behavior; (ii) the linear extrapolation of the 
upper and lower bands appeared to be misaligned , with 
an energy shift corresponding to the gap. 

In order to account for these unconventional features, 
an alternative scenario was proposed in Ref. [24j in terms 
of a phenomenological “massless” gap model with band 
dispersion Ek.± = Ed ±[A/2 + hvk]. The anomalous fea¬ 
tures (i)-(ii) could be naturally reproduced assuming a 
self-energy of the form Ek = A [a x cos 9+ a y sin 8 ], where 
9 = arctan(/cj / /fc a; ). This phenomenological model was 
also shown to provide precise predictions (e.g. on the 
density of states) which could be experimentally checked 


The aim of the present paper is twofold. On one hand 
we would like to provide experimental evidence of the va¬ 
lidity of the massless gap model in CVD grown graphene. 
On the other hand we discuss in details at the theoret¬ 
ical level the possibility of the appearing of a massless 
gap as result of a tendency towards a second-order phase 
transition. We believe that a full understanding and con¬ 
trolling of the unconventional properties of the massless 
gap in graphene on substrates can open new perspectives 
in the bandgap engineering in graphene-based materials. 


II. MASSIVE VS. MASSLESS MODEL 

Before addressing a quantitative analysis of the experi¬ 
mental ARPES dispersion, and a deeper discussion about 
the possible origin of the anomalous bandgap features 
there observed, we summarize briefly here the compari¬ 
son between the massive and massless gap models. 

We first write the non-interacting Hamiltonian in the 
form: 

K = Ed + Hvk ( jj, 6 ) , (2) 

where k are the momenta relative to the K point. The 
eigenvalues of Eq. ©, E± = Ed ± hvk , describe the 
Dirac cone dispersion, while the non-trivial dependence 
of © on the angle 9 accounts for the chiral properties of 
the eigenstates. 

A massive gap is induced if the A and B sublattices are 
electrostatically inequivalent, for instance as an effect of 
the substrate. Such inequivalence is formally taken into 
account by the additional term 

H A = Ed + =E D + ( ^ _^ /2 ) , (3) 


as in Eq. ©, resulting in the well-known dispersion 


E£ s ± = E d ± \J(hvk) 2 + A 2 /4. (4) 


Two regimes can be identified in this energy disper¬ 
sion. At high-energies (k > A/hv) the effect of the 
gap is negligible and the band dispersion recovers the 
normal linear behavior Ej~.± ~ Ed ± hvk. At low- 
energies (k < A/hv) the band acquires a parabolic shape 
Ek,± ~ Ed ± [A/2 + h 2 k 2 /2m e ff \, with m e s = A 2 /2v. 
It is worth noting that in this regime the chiral structure 
of the band dispersion is strongly affected by the opening 
of the gapH. 

The onset of a massive gap in the low-energy elec¬ 
tronic structure of graphene has been discussed not only 
as due to the interaction with the substrate, but also 
as a possible effect of a spontaneous excitonic instability 
induced by the long-range Coulomb interaction [26l - 14l| . 
In this latter context, the term © should be regarded 
as an order parameter , rather than an intrinsic property 
of the system. For real graphene systems, however, the 
dimensionless Coulomb coupling constant results to be 
much smaller than the critical value required for the ex¬ 
citonic instability. On the other hand, in the absence of 
such excitonic broken symmetry, in the normal state the 
long-range Coulomb interaction has been shown to be 
accounted , at the leading order, by a self-energy term 


Ek.Cou 0^ 


( 0 
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k ln(IC/k)e~ iB \ 
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( 5 ) 


Note here the much weaker dependence on k , than the 
bare Dirac Hamiltonian. Indeed, the conventional linear 
behavior in k is here counterbalanced by a logarithmic 
divergence ~ ln(fc) induced by the long-range Coulomb 
interaction. Note finally that in this case the chiral struc¬ 
ture is preserved down to k —> 0. 

The massless gap model investigated in this paper, al¬ 
ternative to the massive model described by Eq. ©. 
resembles features from both © and ©, and it can be 
mathematically described by a self-energy term [24] : 


E m l OC 
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where the off-diagonal terms do not scale to zero for k —> 
0, but they still preserve the full chiral properties. As 
mentioned above, the characteristic energy-momentum 
dispersion results to be 

E^ ± = Ed ±[A/2 +hvk]. (7) 

Note that this models accounts in a simple way for 
the anomalous features described in the introduction, 
namely: i) the Dirac cones appear to be just split by 
the gap A without affecting their conical shape; ii) the 
perfectly linear behavior of the upper (lower) band does 
not extrapolate onto the linear behavior of the lower (up¬ 
per) band. 
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As we are going to see, such anomalies can permit to 
identify in a clear way one gap model rather than the 
other one in the experimental data. 


III. MASSLESS GAPPED DIRAC MODEL: AN 
ARPES EVIDENCE 


In this section we present a detailed analysis of some 
recent ARPES measurements on single-layer graphene on 
Ir(lll) decorated by Ir clusters and Na alkali metals. We 
show how these measurements display, in correspondence 
with the opening of a large energy gap 0 , an anomalous 
behavior of the Dirac cones compatible with the massless 
gap model. 

The experiment was performed at the VUV- 
Photoemission beamline on the Elettra storage ring in 
Trieste under ultrahigh vacuum conditions. The Ir(lll) 
crystal was cleaned by cycles of Ar + ion sputtering and 
annealing at T = 1500 K. Surface order and cleanli¬ 
ness of the sample were checked by low-energy electron 
diffraction and core level photoemission measurements. 
Graphene was grown by thermal cracking of ethylene 
(C 2 H 4 ) on the Ir substrate held at 1300 K for an exposure 
of 100 L (1L corresponds to an exposure of 10 ~ 6 rnbar 
for Is). Under these experimental conditions, graphene 
had lattice vectors aligned in plane to those of the sub¬ 
strate and displayed a Moire pattern originating from 
the interface lattice mismatch [44j |. Tiny Ir amounts were 
evaporated from a current heated thin film plate (0.5 mm 
width, 0.1 mm thickness) at an evaporation rate of about 
2.0 x 10 -4 ML/s, as determined by core level photoemis¬ 
sion measurements 00 . Ir evaporation was performed 
at a temperature of 350 K. The deposition of Ir results 
in the nucleation of size selected Ir clusters on the hep 
regions of the Moire supercells (see Fig. [T|for a schematic 
sketch) 0 ]. The process was finely controlled to satu¬ 
rate the hep regions of the Moire-derived graphene su¬ 
percells with one Ir cluster and to avoid cluster percola¬ 
tion. The optimal Ir coverage was determined to be 0.15 
ML (monolayer) Ir, with reference to the density of an 
Ir(lll) plane fl7t44 |. Additional Na was evaporated from 
a commercial getter source at room temperature and ad¬ 
sorbed to fill the residual exposed graphene surface. The 
resulting system (Na/Ir/G) was characterized by means 
of angle-resolved photoemission, as reported in Ref. 0 . 
The valence band photoemission analysis was carried out 
by a Scienta R4000 electron energy analyzer with 120 eV 



• Ir C Na 

FIG. 1: Schematic sketch of the Na/Ir/G sample. 
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FIG. 2: (a) ARPES intensity map of the electronic disper¬ 
sion of Na/Ir/G along the direction prK. The light green 
points represent the electronic dispersion determined by the 
momentum distribution curves, while white points the elec¬ 
tronic dispersion obtained by energy distribution curves. The 
vertical white double arrow marks the magnitude A of the 
gap. The location of the cut pTK within the Brillouin zone is 
shown in the inset of panel (b). (b) Momentum dependence 
of the average l?M,k (see text). 


photons, close to the Cooper minimum of the Ir5d lev¬ 
els. The spectra were collected at 120 K with energy and 
angular resolution of 30 meV and 0.3°, respectively. 

The ARPES intensity map, along the direction per¬ 
pendicular to TK passing through the K point (pTK) 
is shown in Fig. [2Ji, where the white and green sym¬ 
bols represent the position of the maxima in the energy- 
distribution-curves (EDCs) and in the momentum- 
distribution-curves (MDCs), respectively. The Fermi mo¬ 
mentum along the pTK direction is kp ~ 0.19 A -1 that, 
taking into account the trigonal warping of the bands, 
gives an electron doping of n « 0.028 electrons per C 
atom. From the EDC at k = 0 we estimate a bandgap 
A « 0.8 eV. 

Two remarkable features stand out from the ARPES 
data: (*) both conduction and valence bands present a 
linear dispersion throughout the probed k-space region, 
except within a very small region k < 0.026 A -1 . This 
range is thus much smaller than that predicted by the 
massive gap model (k < A/hv ~ 0.075 A -1 ); ( ii ) the 
linear extrapolations of the upper and lower branched 
do not match, at variance with the expectations of the 
massive gap model. We note additionally the lack of 
any diamond-like features close to the K point, which 
would be typical trademarks of plasmaronic coupling 0- 
[23i | . On the other hand, plasmaron signatures in ARPES 
data have been so far reported only for epitaxial graphene 
on SiC, whereas they are absent in graphene grown on 
metallic substrates, most probably as a consequence of 
screening effects. 

Such band structure of the Na/Ir/G system appears to 
be incompatible with the standard massive gap model, 
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FIG. 3: Comparison between the experimental dispersion for 
Na/Ir/G and the model predictions. Black filled diamonds 
and red empty circles represent EDC and MDC data, respec¬ 
tively. Solid blue lines the predictions of the massless model 
with A = 0.3 eV, and blue dashed lines the corresponding ex¬ 
trapolations of the Dirac cones. The dotted green lines show 
the massive gap model with A = 0.4 eV. In both theoretical 
cases we have used Ed = —1.3 eV and hv « 5.3 eV/A. 


whereas they could be naturally reproduced within the 
context of the massless gap model. 

In order to address at a more quantitative level this 
issue, we provide here a detailed careful analysis of the 
ARPES data for Na/Ir/G. Before testing the different 
models, it is necessary to estimate the basic parameters 
which are independent of the model itself, e.g. the Fermi 
velocity v and the the Dirac energy E d in the absence of 
the gap. 

We estimate Hv from the linear behavior of the lower 
band in the energy window from —3.5 eV to —1.7 eV, 
which yields Hv ~ 5.3 eV/A. This value is consistent with 
the speed value derived from the upper band. We found 
this value was also perfectly compatible with the linear 
slope of the upper band. A rough estimate Ed = — 1.3 
eV of the Dirac energy in the absence of the gap is ob¬ 
tained as the midgap point at k = 0 (E k= o,+ = —0.9 eV, 
Ek=o- = —1-7 eV). This estimate can be cross-checked 
with the average energy between the upper and lower 
band E M ,k = (E k = o,_ + E k= o,-)/2 for generic k, dis¬ 
played in Fig. [2jo. The result shows a negligible mo¬ 
mentum dependence, resulting in i?M,k = —1.3 eV. Note 
that the negligible dependence of Eyi,k on k signalizes 
an almost perfect symmetry of the experimental upper 
and lower bands with respect to the Dirac energy. This 
is compatible with both the massless and massive gap 
models, but it rules out plasmaronic effects which would 
result in an asymmetry and in a finite k dependence of 

The unbiased determination of Ed and v allows to test 
the validity of the different models on the experimental 
data. Fig. [3] compares the band dispersion extracted 


from the momentum and energy distribution curves with 
the prediction of massless and massive modes. The best 
agreement is found for the massless model with A = 0.3 
eV. Note the mismatch between the upper and lower 
Dirac cones, that reproduce perfectly the experimental 
data. The same value of the gap can be indeed also in¬ 
ferred from such misalignment. On the other hand, no 
satisfactory fit was possible for the massive model. We 
show here thus the predictions of the massive model for 
A = 0.4 eV, that reproduces the value at K of the upper 
and lower bands. 


IV. MICROSCOPIC MODELS 

In the previous section we have shown how a care¬ 
ful analysis of the band anomalies observed by ARPES 
in in Na/Ir/G points out towards a phenomenology of 
the graphene bands compatible more with massless gap 
model rather than with the standard massive one. It is 
clear that the different models represent a trademark for 
different underlying physical processes. In this Section 
we address thus this issue, in order to clarify which kind 
of physical processes can be responsible at the micro¬ 
scopic level for the massless Dirac phenomenology, with 
the final aim to control and engineering such processes. 

Along this perspective, a brief summary of the mi¬ 
croscopic physics related to the massive model can be 
useful. As mentioned in the Introduction, the massive 
Dirac model essentially stems from the presence of a lo¬ 
cal potential term oc a z which makes the two carbon 
atom sublattices inequivalent. This potential can be nat¬ 
urally induced by the interaction of graphene with the 
substrate. The term A a z in the Hamiltonian appears 
then as an external field arising from the environment. 
Alternatively, the onset of a mass term oc d z has been 
discussed within the context of a spontaneous symme¬ 
try breaking, where the driving mechanism is the un¬ 
screened long-range Coulomb interaction [2(il-(4ll]. The 
term A(k)tf z can be viewed thus in this scenario as an 
anomalous self-energy term, i.e. as an order parameter , 
arising from many-body effects. It should be however 
mentioned that in this case the mass term A(k) appears 
to be momentum dependent, i.e. divergent for k —> 0 
rather than k-independent. 

From a general point of view, since the presence of a 
linear Dirac cone in ideal flat graphene is dictated by 
symmetry reasons, it is clear that the opening of a gap 
and the splitting of the upper and lower bands, both in 
the massive and in the massless model, must be associ¬ 
ated with some kind of breaking symmetry. We have in¬ 
deed just seen how the symmetry breaking (spontaneous 
or not) responsible for the massive model is the equiva¬ 
lence between the local potentials of the two sublattices. 
In Ref. [13 the general constraints that a phenomeno¬ 
logical self-energy must satisfied to give rise to a massless 
Dirac model have been discussed. It is on the other hand 
intriguing to discuss in the present work the conditions 
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that could make possible the appearance of a self-energy 
m as result of a tendency towards a spontaneous sym¬ 
metry breaking. 

It is clear that the effective breaking of symmetry, 
in real materials, can be prevented by many causes. 
For instance, as we will discuss below, a good candi¬ 
date for a quasi-breaking of symmetry is the long-range 
Coulomb interaction. In real materials, and in partic¬ 
ular in Na/Ir/graphene, such interaction is expected to 
be screened resulting in a cut-off of the singular behav¬ 
ior for small exchange momenta q. Asymptotically, this 
should correspond to “regularization” of the interaction, 
and thus to a linear behavior of the dispersion (although 
with a very steep slope) for momenta k close to the K 
point. For larger k, however, the tendency towards an in¬ 
stability can still be reflected in a quasi-gapped dispersion 
which can be compatible with the experimental observa¬ 
tions. From this point of view, we present the analysis of 
the onset of a spontaneous symmetry breaking just as an 
illustrative example for the gross features of the result¬ 
ing phenomenology. In similar way, we do not want to 
provide here a quantitative evaluation of the possible crit¬ 
ical couplings required to induce a massless gapped case, 
rather to compare, at a qualitative level, this strength 
of this instability with the strength of the corresponding 
instability associated with the spontaneous generation of 
a mass according the model ©• In order to have thus a 
direct and clear way to compare the two possible insta¬ 
bilities, we address here this issue within the context of 
a mean-field solution, keeping in mind that a more com¬ 
pelling high-order analysis is needed for a quantitative 
investigation. 

In order to address this issue in the most convenient 
way, we introduce the Nambu notation where the Hamil¬ 
tonian of ideal free-standing graphene can be written as: 

H = ^k,CT-^o( k )V’k,er) (8) 

k ,cr 

where a is a global index taking into account both spin 
and valley (degenerate) degrees of freedom N s = 2, 
N v = 2), ipk a spinor in the sublattice basis 

= ( c k,<x,A> c k,a, B ). and where 

Ho(k) = Hvk [a x cos 9 + d y sin 9] . (9) 

We assume a long-range Coulomb interaction which 
can we written in the momentum space as: 

= a°) 

c t,ct' k,p,q 

where p q = Ek,cr V’k+q.crV’k.a, and where V(q) = 
27re 2 [l — <5(q)]/eo«|q|. The term oc [1 — <5(q)] in V(q) 
takes into account the subtraction of the positive charged 
background, while k is the dimensionless relative in-plane 
dielectric constant. 

The phenomenology of the normal state self-energy 
arising from such many-body interaction, as well as of 


several possible anomalous self-energies associated with 
different spontaneous symmetry breaking, has been dis¬ 
cussed in an extended way using different accurate tech¬ 
niques, as for instance renormalization group (RG) or 
static random-phase approximation (RPA). In general, 
in order to have an accurate estimate of the magnitude 
of the corresponding self-energies (and of the minimum 
strength of the interaction in the case of spontaneous 
symmetry breaking), a self-consistent screening of the 
quasi-particle excitations and of the the effective screened 
Coulomb interaction is required. 

On the other hand, a qualitative insight on the ma- 
trix/momentum structure of the possible anomalous self¬ 
energies is already possible employing a simple mean-field 
treatment. 

Neglecting the Hartree term, which does not play any 
role in the present context, and focusing on the exchange 
(Fock) term, we can write the generic self-energy as: 


E(k) = TV- . V(k P) . -. 

f^iu n I-H 0 ( k)-E(k) 


( 11 ) 


Note that, since the interaction V(q) is peaked at small 
q = k — p, the intra-valley scattering is dominant, sup¬ 
porting the validity of the low-energy model versus a full 
tight-binding treatment. Note also that at this level of 
approximation (unscreened static Coulomb interaction) 
the self-energy does not depend on the frequency iui n . 
Within this scheme is thus convenient to define an effec¬ 
tive mean-field hamiltonian H{ k) = Ho(k) + E(k), and 
to write a recursive equation: 

H(k) = H 0 (k) + Tj2 . (I 2 ) 

p,„ - H(k) 


Since R(k — p) depends on the angular variables 9 &, 
9 p essentially through only cos (9k — 8 P ), one can see that 
H{ k) f3(k)[a x cos9 + d y sin61 Apart numerical pref¬ 
actor, the momentum dependence of /3(fc) on k can be 
caught already by a perturbative analysis. Replacing 
H( k) with Ho(k) in the right side term of Eq. (TT21) . 
or, equivalently, neglecting E(k) in the right side term of 
Eq. (ED, we obtain: 


/3(k) = Hvk 


1 + s 7(t) 


(13) 


where a = e 2 /eo hv is the coupling constant for suspended 
graphene, K, is a high-momentum cut-off limiting the va¬ 
lidity of the Dirac model, approximately determined by 
the size of the Brillouin zone, and where 


7(fc) 



d9 x cos 9 
27T y/\ + x 2 — 2a; cos 9 


(14) 


A careful analysis of (fHl) shows that at the leading or¬ 
der 7 (fc) ss ( 1 / 2 ) ln(/C/fc)+const., a well known result. 
Note that, within this context, the self-energy associ¬ 
ated with the massless gap model is formally described 
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by (3(k) =const., i.e. by j(k) = 1/k. It is clear that such 
analytical dependence is not spontaneously generated by 
the self-consistent solution of m in the normal state, 
and it must be regarded as an order parameter of a pos¬ 
sible second order phase transition. In order to address 
this possibility, we rewrite Eq. m as: 

H(k)-W[k,{H(p)}} = H 0 (k), (15) 

where IU[k, {A}] is the functional W[k, {X}] = 

n U(k — p )/[ioo n I — X). Broken symmetry phases 
can be generated when the equation 

0(k) - W[k, {H(p) + <£(p)}] = 0 (16) 


It is remarkably relevant to notice that, contrary to 
Eq. (1141) . the integral over p is well behaved for both 
k —> 0 and /C —> oo. In this limits we can thus properly 
look for the self-consistent solution (3 m i(k) = A m i. We 
obtain thus the susceptibility equation: 

ry 

1 = TTX ml (19) 


where 


Xml 



dd cos 9 

27r y/1 + x 2 — 2xcos 9 


2- V3/2 


1.13. 


( 20 ) 


admits a non-trivial solution ^(k) ^ 0. We are here inter¬ 
ested only in investigating the instability of the normal 
state towards a broken symmetry phase. Within this 
context we can expand thus the right side term of Eq. 
m at the linear order with respect to <^(k). At T = 0 
we get thus the susceptibility equations: 

m = (i7> 


where E( p) = det[I0(p)] is the energy dispersion in the 
normal state. Within the spirit of a mean-field approach 
we can neglect in E(p) the self-energy many-body effects 
driven by the Coulomb interactions, and approximate 
E(p) ~ det[.ffo(p)] = fiv |p|. Eq. El) can be used to 
investigate the instability towards both a massive as well 
as massless gap. In the first case </>(k) oc d z , which clearly 
breaks the symmetry of the original Hamiltonian, while 
in the second case </>(k) oc a x ,a y . Since terms oc d xi a y 
are already present in the normal state, it is not straight¬ 
forwardly apparent at this level the nature of the symme¬ 
try breaking. On this regards it should be noticed that in 
the massless gap model the symmetry breaking is more 
specifically associated with the anomalous k-dependence 
of the order parameter, in similarity with q = 0 flux 
phases discussed in different contexts, rather than with 
the Pauli structure. 

To make clearer this point, we employ the massless 
gap model 0 m i(k) = /\ m \[a x cos 9 + d y sin0] as an ansatz 
for the right side of Eq. (1171) . It is easy to check that 
the equations for cr Xl a y result to be decoupled and de¬ 
generate. From a careful analysis of the angular vari¬ 
able, it is also straightforward to check that the result¬ 
ing (^(k) in the left side of Eq. (fT71) appears to be in¬ 
deed 4>(k) = (3 m \(k)[d x cos9 + <7 y sin0], where the pref¬ 
actor ,0 m i (k) will be determined by the effective self- 
consistency of Eq. m- Using the explicit expression 
for the Coulomb interaction, we obtain thus: 


/3 m i(fc) = A 


a 

ml 2 k 


x 



dp 


[ cos0 (is) 

Jo \Jk 2 + p 1 — 2 kp cos 9 


For free-standing graphene with k = 1 we can predict 
thus a finite critical coupling a™ s ~ 1.76 above which the 
system undergoes spontaneously a transition towards a 
massless gap phase. 

We should stress that such value of a™ s has been de¬ 
rived within the lowest order approximation, and a more 
compelling analysis is needed for a reliable quantitative 
estimate. It is however instructive to compare such value 
with the critical coupling required to induce, at the same 
level of approximation , a spontaneous generation of mass. 
Such instability can be determined by using the ansatz 
</> ms (k) = [A ms /Vk\a z in Eq. (El). After few straight¬ 
forward steps, we obtain the susceptibility equation 

a 

1 = 7T“Xms (21) 


where 


Xms 


r°° dx_ 

Jo 

r 4 (i/4) 

47T 2 



d9 

2n 


\/l + x 2 — 2x cos 9 


~ 4.38, 


( 22 ) 


corresponding, for k = 1, to the well-known result aU s ~ 
0.46 glial. 

Although the critical coupling a™ 1 for the spontaneous 
generation of massless gap results to be larger than the 
critical coupling for the massive gap, it remains still of 
the same order of magnitude of a™ s within the same level 
of approximation. We know however that a™ s is signif¬ 
icantly increased when higher order renormalization ef¬ 
fects (dynamical screening, self-energy effects, ...) are 
taken into account, with possible a™ s as large as a™ s ~ 2, 
of the same order thus of the present estimate of a™ 1 . Be¬ 
comes thus crucial, in order to assess the dominant insta¬ 
bility of the normal state, to include high-order effect also 
for the massless gap case, to permit a quantitative com¬ 
parison. In addition, it should not be excluded a possible 
interplay, in the broken symmetry phase, between the two 
order parameters, where one could favor and make more 
stable the other one. This scenario is suggested by the 
simple analysis of the experimental data of Fig. [3l where 
the massless gap model reproduce well the experimental 
dispersion in a broad energy range, but where a finite 
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parabolic behavior is nevertheless present at very small 
k. We have already discuss how such weak parabolic be¬ 
havior cannot account in any case for the large observed 
gap. The experimental dispersion could be thus compati¬ 
ble with the possible coexistence of the both qualitatively 
different gaps, with the massless (larger) one determin¬ 
ing the high-energy range, and the massive (smaller) one 
accounting for the residual parabolic behavior at small 
k. Further analysis in this direction is however needed to 
investigate this issue. 

In conclusion, in this work we have provided, using 
angle-resolved photoemission spectroscopy, a compelling 
evidence of the presence of an unconventional gap in 
CVD grown graphene, properly described by a mass¬ 
less gap model. We have further investigated, on a 
theoretical ground, the microscopic conditions that can 
give rise to a spontaneous phase transition towards such 
massless gapped phase. We have shown that, for sus¬ 
pended graphene, a critical value of the Coulomb cou¬ 
pling is needed, with larger but comparable values with 


the ones required to induce a breaking symmetry towards 
a massive map. This analysis, stricty valid for suspended 
graphene, can suggest a possible path to induce such un¬ 
conventional gap also in CVD graphene, although in this 
case specific features of these materials (screening of the 
Coulomb interaction, interference with the Moire pat¬ 
tern) must be taken into account. The fully understand¬ 
ing and controlling of the unconventional properties of 
the massless gap in graphene can open new perspectives 
in the bandgap engineering in graphene-based materials. 
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